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BETHE-PEIERLS-WEISS APPROXIMATION AND A MODEL 
FOR FERROMAGNETIC THIN FILMS 
by Gale Fair 
Lewis Research Center 

SUMMARY 

The cluster approximation developed by P. R. Weiss to treat the problem of the 
Heisenberg ferromagnet is used to calculate the ferromagnetic properties of a two- layer 
rectangular lattice, a simple cubic lattice, and a body-centered cubic lattice, for spin 
1/2 atoms. The two-layer rectangular lattice is considered as a model for a simple cu- 
bic ferromagnetic thin film, and the results are compared with the properties of a com- 
plete simple cubic lattice. The results obtained indicate that such a two-layer lattice 
shows ferromagnetic behavior, with a Curie temperature about two-thirds that of the sim- 
ple cubic lattice. The magnetization, internal field, energy, specific heat, and paramag- 
netic susceptibility are calculated and compared for the three lattices. One feature that 
appears in this calculation for the two-layer lattice, as previously obtained for the sim- 
ple cubic case, is the occurrence of an M anti-Curie M temperature; a second, lower, 
temperature at which the spontaneous magnetization vanishes. 


INTRODUCTION 

The existence of ferromagnetism for some metals (e. g. , iron, nickel, cobalt, and 
some of their alloys) is a phenomenon for which no completely adequate theoretical ex- 
planation has yet been given. The actual interaction mechanism that gives rise to ferro- 
magnetism is not fully understood. In practice a Hamiltonian is assumed for a model of 
this interaction, such as the Heisenberg Hamiltonian, which will be used in this work. 
Once the Hamiltonian is specified, however, the problem is not finished. For example, 
in the case of the Heisenberg model, the Hamiltonian eigenvalue problem cannot be 
solved exactly for a crystal. To determine the thermodynamic properties of a system 
from this model, several techniques have been developed for making approximate calcu- 
lations. 



One of these approximation methods is the method of spin waves first introduced by 
Bloch (ref. 1). The original Bloch theory predicts the existence of a Curie temperature 
for a three-dimensional lattice but finds no such Curie temperature for a planar, two- 
dimensional lattice. This method has also been applied to the case of thin ferromagnetic 
films (refs. 2 and 3), to determine the effect of film thickness upon the spontaneous mag- 
netization for films in which the thickness of the lattice is quite small in comparison with 
the planar dimensions. The results of spin-wave theory applied to finite lattices (refs. 2 
and 3) indicate that such systems, at sufficiently low temperatures, show ferromagnetic 
behavior, even to the limit of lattices that are only one atomic layer in thickness. 

The primary motivation of the present analysis was to determine whether these re- 
sults of spin-wave theory were obtainable with other methods of dealing with ferromag- 
netic phenomena. The method selected was the Bethe-Peierls-Weiss approximation (or 
BPW method), which is basically a high-temperature approximation for treating the 
Heisenberg Hamiltonian. This method is applied herein to a new model for a ferromag- 
netic thin film. This thin- film model consists of two plane quadratic (square) lattices 
parallel in space with the distance between the planes equal to the interatomic distance 
within each lattice. This means that each atom will have five nearest neighbors. 

In addition to the analysis of this thin-film model, a more complete investigation is 
made than was available of the properties of other simple lattices. This permits com- 
parison of the results for the thin-film model, not only with the low temperature spin- 
wave calculations, but also with the properties of bulk samples treated with the BPW 
method. 

To clarify the nature and source of the BPW method to be used, it is necessary to 
review briefly the more elementary Ising model, with its approximate methods of solu- 
tion, as well as the Heisenberg model. 


Ising Model 

The Ising model represents a method devised by Ising (ref. 4) (and subsequently de- 
veloped further), to treat a class of cooperative phenomena. This model has been used 
in the theory of ferromagnetism, the "lattice gas, " order- disorder in binary alloys, and 
other problems involving interactions between atoms situated upon lattice sites. It is as- 
sumed that there is an interaction energy between a lattice atom and each of its nearest 
neighbors in the lattice. For example, in ferromagnetism this energy is to depend upon 
the relative orientation of the spins of the atoms; for order-disorder it depends upon 
whether the neighboring atoms are like or unlike species; and so forth. The problem is 
to use statistical mechanical techniques to calculate the thermodynamic properties of the 
system from a knowledge of its total energy. 
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As a model for ferromagnetism, the Ising model is too highly idealized to predict 
quantitative results for a real ferromagnet. At most, this model may be used as a guide 
for some of the qualitative aspects of a ferromagnetic system. In this vein, Onsager 
(ref. 5) has given a rigorous calculation for the two-dimensional Ising system and finds 
that such a lattice has a Curie temperature. 

The total energy for the Ising model is determined by the Hamiltonian for the system. 
This Ising Hamiltonian has the following form: 

“is -- 2 Z>i s r B Z s i <») 

i>3 i 

The first summation represents a double sum, one over all spins S^, and the second over 

all the nearest neighbors to S., denoted by S-, where S, represents, for ferromagne- 

1 1 1 th 

tism, the value of the scalar spin of the atom that is on the i lattice site, and e is the 
interaction energy between S. and Sj. The second term represents the interaction be- 
tween the spin system and an external magnetic field proportional to B. 

Among the approximate methods used to solve this problem, the one of particular in- 
terest here was developed by Bethe (ref. 6) for the solution of the order -disorder prob- 
lem in binary alloys. It was later shown by Peierls (ref. 7) that this method also could 
be extended to the treatment of ferromagnetism. Basically, this approximation consists 
of treating rigorously the interaction of a single atom with its nearest neighbors. The 
remainder of the system is then replaced by an equivalent internal field that interacts 
only with this first shell of nearest neighbor atoms. This problem can be solved exactly 
and the thermodynamic properties determined from the corresponding partition function, 
which includes the internal field. Since there must be no difference between the proper- 
ties of the atom chosen as central and those for an atom in the first shell, a self- 
consistency condition is used to eliminate the internal field. In practice, this is done by 
requiring that the average magnetic moment of the central atom in terms of the external 
field and the average magnetic moment of an atom in the first shell in terms of the inter- 
nal field be calculated and set equal. This calculation explicitly determines the internal 
field. 

This procedure may be extended to a series of approximations. At each step one 
more shell of atoms is included in the rigorous calculation, the remainder still being re- 
placed by the internal field. This procedure is extremely difficult beyond the first step, 
however. 

The method used in this report is an extension of the Bethe-Peierls approximation 
developed for the Ising model (as described previously), to the Heisenberg exchange 
Hamiltonian. 
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Heisenberg or Exchange Interaction Model 

The Heisenberg model (ref. 8) for ferromagnetism was developed as a consequence 
of the exchange interaction of quantum mechanics. This exchange interaction between a 
pair of spins arises from the requirements of the Pauli exclusion principle and the 
Coulomb interaction and may be thought of as a result of the overlapping of the wave func- 
tions describing the atoms. 

Similar in form to the Ising Hamiltonian, the Heisenberg Hamiltonian for a system 
of spins that occupy lattice sites is given by 

*H-- ! I j vv B E s k (2 > 

i>j i 

plus an unimportant constant term involving only J. The summation is over nearest 
neighbor pairs, J is the exchange interaction coupling constant, and S- is the z- 

J.'U 

component of the 1 spin. The vector spins here are to be treated as quantum- 
mechanical operators. 

Among the methods used to treat equation (2) is an extension due to P. R. Weiss 
(ref. 9) of the Bethe-Peierls approximation for the Ising model. This method (referred 
to as the Bethe-Peierls-Weiss approximation, or simply the BPW method) is the one 
used herein. The details of the calculation will be given in the next section of this report, 
but it is basically the Bethe-Peierls method applied to the exchange interaction. One 
atom is chosen as the central atom, and its interaction with the first shell of nearest 
neighbor atoms is calculated exactly; the rest of the system is replaced by an internal 
field that acts only upon this first shell. The corresponding energy eigenvalue problem 
can be solved and the resulting energy levels used to calculate the partition function for 
the system. From this partition function and the self-consistency condition on the aver- 
age magnetic moments of the atoms, the ferromagnetic properties for the lattice may be 
determined. 


THEORY 

Energy Levels and Partition Function 

The cluster Hamiltonian for the BPW approximation is obtained from the Heisenberg 
exchange Hamiltonian (eq. (2)), which is for the entire crystal, by the procedure sketched 
in the introduction. The exchange interaction between the spin of the chosen central atom 
Sq and the resultant spin of the first shell of nearest neighbors Sj is extracted from the 
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exchange term in equation (2). The remaining interactions of the first- shell atoms with 
their nearest neighbors and with the applied magnetic field are then written as an inter- 
action with an effective internal field Bj. The external applied field Bq defines the 
z-axis of the system. In this approximation the internal field Bj, as described before, 
is assumed to be in the same direction. Aside from the exchange interactions, the spin 
of the central atom interacts only with the externally applied field, and only the first 
shell of atoms interacts with the internal field. With these conditions, the Hamiltonian 
for a cluster in which there are no interactions between atoms within the first shell has 
the following form: 


3C = -2JS q • Sj - gi3S 0 z B 0 - g^S lz B 1 


( 3 ) 


where 

J exchange integral 

Sq z z- component of spin of central atom 

Sj z z-component of resultant spin of first shell 

Bq magnitude of external field 

Bj magnitude of internal field 

g Lande spin factor 

18 Bohr magneton 

For convenience in notation, g/3Bg will be replaced by Hg and g/3B^ by Hj. 

The eigenvalues of Sq • Sj may be calculated from the vector model of angular mo- 
mentum as follows: 


S = Sq + Sj 


Therefore, 


S • S — Sq • Sq + Sj • Sj + 2 Sq • Sj 

thus, 

2Sq • S x = S- S- Sq. Sq - S 1 - Sj (4) 

Now define the quantity E„ v as the expectation value of the exchange interaction, as 
follows: 
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E ex' - 2J ( S 0' S l> 


( 5 ) 


Thus E^ will have the values E_ and E + for S = Sj - ^ and S = Sj + ^, respectively. 
Then, from equation (4) 


E_ = J(Sj +1) S = Sj - i 


E + =- JS i S=S l + f 


1 


( 6 ) 




when the spin per atom is 1/2, the only case considered here. 

The representation used is one in which the exchange energy is diagonal. Since 
Sq z + Sj z = S z = m, the z- component of the total spin, by rearrangement of the field- 
dependent terms in equation (3), 

®0z E 0 + ®lz H l - S 0z H 0 + < m - S Oz> H l = “Hi - S 0 z< H l - Hq) 

In this representation, the energy levels when m = ^ are 

W = " JS 1 T ~ ( 2s i H i + h q) 

When | m | ^ the matrix for the Hamiltonian is (see appendix) 


a 




B 

V 


where 


A a = J (Sl + 1) - m Hl - -BL- ( Hl - H 0 ; 


A p = -JS 1 -mH 1+ -^-(H 1 -H 0 ) ^ 


B = — 


4rn 


11/2 


( 7 ) 


(2S 1 + 1)‘ 


(Hj - Hq ) 
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For convenience, the following definitions are made: 


xsH i- H o 


Of 


m 


2Sj + 1 


The determinantal equation for the eigenvalues then becomes 

1 I 

E - mHi - ax - A ± (l - 4 a*) x 

1 2 ' 

,1/2 


— (l - 4d! 2 ) x E + - mHj + ax - X 


= 0 


( 8 ) 


The resulting quadratic equation for A^., the eigenvalues, may be solved exactly and 
yields 


E + E 


E - E 


k 2 mHl * + 2 ~ & + |X(4a + ^ 


X,. = 


1/2 


( 9 ) 


where 


E + - E_ <1(2^ + 1) 

The term £x is proportional to the ratio (H^ - Hg)/J. Ultimately H Q will be set equal 
to zero; in the neighborhood of the Curie point, the internal field is small, so that 
the square root may be expanded in terms of the ratio £x. After expanding to fourth 
order and rearranging terms, the energy levels are 

x k* W k= -Wl +c k< H l - H 0> +d k< H l - H 0> 2 + e k< H l - H 0> 3 + V H 1 - H 0>“ *-1. 2 

( 10 ) 


where 


■^1,2 E =f 

b i = b 2 = -m 
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C 1 , 2 = T 


m 


2Sj + 1 


di o _ ± 

A ’^ 4J(2S^ +1) 


1 - 


4m" 


(2^ + 1)' 


e i o 

- t 2 


m 


2J (2Sj + 1)' 


4m 


(2S^ + 1)' 


f l,2 = * 


16J 3 (2S 1 + l) 3 


, 24m 2 80m 4 

1 + 


(2Sj_ + l) 2 (2^ + l) 4 


Terms to fourth order in the field strengths are kept for the specific heat calculation 
although only the second-order terms are needed in explicit form for the determination 
of the Curie temperature. 

The partition function Z for the cluster, written in terms of the energy levels W^, 
then takes on the following form: 


-z<^>z 


-W 1 (S 1 ,m)/iT 
e A + 


m 


I 


m . 


-W 2 (Spin)/iT 


( 11 ) 


The first summation is over the value of the total spin of the first shell Sj from Sj = 0 
to Sj = n, where 2n is the number of nearest neighbor atoms within the first shell. The 
summation indicated by m_ is over m from m = to m = and the 

m + sum is from m = - (Sj + to m = +^Sj + The factor u>(2n, Sj) represents the 
degeneracy of each level. This is essentially the number of ways that the first shell 
of atoms may have a total spin of S^. To calculate this factor, let L(2n, S^) stand for 
the number of ways of placing 2n spins, n + up and n - Sj down, upon the 2n lat- 
tice sites, giving a z-component of the total spin equal to Sj, with the total spin of Sj 
or greater. Then 


L(2n, S 1 ) = 


(2n)! 


(n + S^.'C n - Sj)! 


( 12 ) 
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The corresponding factor for the arrangement of n + + 1 up spins and n - Sj - 1 

down spins upon the 2n sites, resulting in a total spin of + 1 or greater, is 
L(2n, Sj + 1) 


L(2n, S x + 1) 


(2n)! 

(n + + l)!(n - S 1 ~ 1)* 


( 13 ) 


The difference between these two factors is just the number of arrangements upon the 2n 
sites that give a total spin exactly equal to Sp with a z- component also equal to Sp 
Thus 


co(2n, S x ) = L(2n, - L(2n, S x + 1) 

(2n); (Ml (14) 

(n + Sj) ! (n - Sj)! (n + Sj + 1)! (n - Sj - 1)! 

Self-Consistency Condition and Curie Temperature 

The expression for the partition function (eq. (11)) still contains the internal field 
explicitly. To eliminate this from the calculation, the consistency condition must be ap- 
plied to the system; that is, the central atom must be identical to any atom in the first 
shell. Explicitly, this means that the average magnetic moment of the central atom must 
be equal to the average magnetic moment of an atom in the first shell. Since the central 
atom does not interact with the internal field, the average magnetic moment is given 

by 


mg = hn 


d In Z 
3B„ 


= g|3^T 


d In Z 


SH 


0 


(15) 


The average magnetic moment of a first-shell atom, which is determined by only the 
internal field, is given py l/2n times the average magnetic moment of the entire first 
shell and is 


m 1 = — <4 t — — = — gj3^T - ln - 
1 2n SBj 2n SHj 


Thus, for the consistency condition, niQ and m 1 must be equal 


(16) 
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2n - _ 

3H 0 3H 1 


( 17 ) 


The solution of this equation for the internal field can be obtained for the case of a 
vanishing external field. Thus, after the derivatives of the partition function are taken, 
Hq is set equal to zero. In evaluating equation (17) odd functions of m must vanish, 
since the sums are taken over an interval symmetric about m = 0. Equation (17) gives 
the following equation determining the internal field explicitly: 

AHj + CHj = 0 (18) 

which for nonvanishing yields 

A + CHj = 0 (19) 

The quantities A and C are defined as 


2 n 


with 


"EE*- 

k=l S ± 

2 n 

'■EE w(2n, 
k=l S 1= 0 


0 


,0 , , u k , °k 


Cj, + — + + 


(kT ) 2 (^T ) 3 


-\/kl 


m 

(1 + Zn)^ [(b k + Cj.) 2 + 2nc k (b k + 

m 

C k s Z ( ‘ 4fk) 

m 

+ 2n)C k = Z 2 [ (i + 2n)d£ + (2 + ^n)\e k + 2(1 + 2n)c k eJ 

m 

2n)c k 3 Z f 2 d k< b k + °k > 2 - 2 nd k< b k + °k> (b k + 2 c k ) ] 

rn 

6(1 + 2n)C 3 s Z K + Ck)4 + 2nCk(l>k + Ck>3 ] 


-(1 + 


( 20 ) 


( 21 ) 
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To evaluate these summations the following identities are used: 


T 

^ 1 = 2T + 1 
m=-T 

T 

Z m 2 = — T(T + 1)(2T + 1) 

3 

m=-T 
T 

m 4 = — T(T + 1)(2T + 1)(3T 2 + 3T - 1) 

15 

m=-T 

It is important to note that in deriving equation (19) one factor of Hj has been can- 
celled from each side. This implies that = 0 is always a solution to the consistency 
condition. This result is analogous to the Bethe-Peierls solution for the Ising problem, 
which also gives such a trivial solution. For a system for which there exists a nonzero 
Curie temperature, this solution will not be of interest. Physically, this zero solution 
corresponds to no long range order in the system. 

The quantities A and C in equation (18), the consistency condition, are functions 
of the temperature, or, more precisely, functions of the ratio x = J/<4t. Thus, Hj is 
determined as a function of x or T. The Curie temperature T c is specified by the 
vanishing of the internal field, or, equivalently, by the vanishing of the quantity A in 
equation (19). Explicitly this means evaluating A, setting the expression equal to zero, 
and solving for x. This equation has the general form 

n 

o>(2n, S<) 4 
1 3 

Sj=0 


IK 





/ i\ / n -(Si+i)*i r s i x -( s i +i ) x ii 

( s i-“)( s i4) (i+n+s i )e J + L e e ir° 


ii 



Equation (23) contains the parameter n, which specifies the type of lattice under con- 
sideration. Thus equation (23) determines the Curie temperature for any lattice that has 
no interactions among the first-shell atoms. For the different values of n correspond- 
ing to the lattices of interest here, the implicit equations for x £ = J/kl c are as follows: 

(1) For the linear chain 2n = 2, equation (23) becomes 

e 3x = 1 - — x (24) 

2 

and the only real, nonnegative solution is x fi = 0. Since this is a solution that is not pos- 
sible physically (infinite temperature), the only possible solution of equation (18) for this 
system is = 0. 

(2) For the quadratic lattice 2n = 4, equation (23) becomes 

l 2 e 2x - 5(8x - l)e x + (32x - 5)e" 2x +— (x - iV 3x = 0 (25) 

25 5 X 5/ 

Again the only nonnegative, real solution is x c = 0, making the only physically possible 
solution of equation (18) = 0. Thus, for this two-dimensional lattice there is no Curie 

temperature in agreement with the Bloch spin-wave theory. 

(3) Similarly, for the simple cubic lattice 2n = 6, equation (23) has the following 
form: 

2 / 3x -4x, 1 / 2x -3xx x/7 2x n -3x\ 1 ,x -2x^ 

— (e -e )+— (e -e ) + — ( — e -9e ) + — (e - e ) 

49 5 7 \2 / 3 

5 / x 1 -2x\ 5 -4x n /0£n 

— x(e — e )+ — xe =0 (26) 

7 V 2 7 7 

One acceptable solution to this equation is x p = 0. 540 or <4 t c = 1. 850 J. 

(4) For the model for a two-layer thin film, a lattice with five nearest neighbors, 
the equation for x fi is 

-1 [fee- 3 * +1 <e 3x - e- 3x )l + i (e 2x - e' 231 ) - * (3e 2x - 5e" 2x ) + * <e x - e' 51 ) - * e x = 0 

36 L 3 J 4 6 4 3 

(27) 

The solution of equation (27) is x fi = 0. 783, corresponding to a Curie temperature given 
by ^T c = 1.278 J. 
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(5) The explicit form of equation (23) for the body-centered cubic lattice (2n = 8) is 
quite complicated, but the numerical solution gives the Curie temperature of 
It c = 2.908 J. 

The previous results show a qualitative agreement with the Bloch spin-wave theory. 
The one- and two-dimensional lattices show no ferromagnetic behavior, that is, no Curie 
temperature. The three-dimensional space lattices have solutions for the Curie temper- 
atures, while the five-nearest-neighbor lattice structure, constructed as an intermediate 
model, also has a Curie temperature, which falls below that for the simple cubic lattice. 


Anti-Curie Temperature 


In addition to the Curie temperature solutions for the consistency condition as dis- 
cussed previously, there also appears a second set of solutions to equation (23) for those 
lattices for which there exists a nonzero solution for x. These second solutions, which 
occur at lower temperatures than T , have been termed anti-Curie points by P. W. 
Anderson (ref. 10), who first noticed their existence. For the lattices for which there 
is no real solution to equation (23), the question does not arise. The simple cubic lattice 
2n = 6 has such an anti-Curie point with the value ^T^ c = 0. 269 J as found by 
Anderson. For the two-layer model 2n = 5 and for the body- centered cubic 2n = 8, the 
second solutions occur at At ^ = 0. 339 J and Ut ^ = 0. 212 J, respectively. 

The existence of a second zero of the internal field, and correspondingly, the vanish- 
ing of the magnetization of the system below T^-, is clearly in disagreement with the 
physical behavior of a ferromagnetic system. 

The source of this difficulty can be traced to the expansion of the energy eigenvalues 
in powers of (H^ - Hq)/J. The region of temperature for which this expansion is expected 
to be valid can be obtained by evaluating the internal field, obtained from equation (19), 
as a function of temperature. For the simple cubic lattice and also for the body-centered 
cubic 2n = 8, there is a singularity in the internal field; the factor C in equation (19) 
vanishes. Numerical evaluation shows that this singularity occurs very near the Curie 
temperature, at a temperature of about 90 percent of the Curie temperature. For the 
two-layer lattice the factor C does not vanish, but it becomes small enough that the in- 
ternal field becomes very large, also quite near the Curie temperature. The appearance 
of the anti-Curie point cannot be taken too seriously since in all cases it occurs outside 
what must be considered as the region of validity of the energy eigenvalue expansion. 

Figures 1 and 2 show the calculated dimensionless internal field and magnetization 
mg very near the critical temperature for the three lattices under consideration. All 
curves show the extremely rapid decrease to zero at the transition temperature, which 
is characteristic of ferromagnetic transitions. 
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Energy and Specific Heat 


Since the specific heat of the system is the temperature derivative of the energy, the 
starting point for this calculation is the specification of the energy per atom of the system 
as a function of the temperature. The exchange term in equation (3) represents the inter- 
action between a central atom and its shell of first nearest neighbor atoms. In a system 
with N atoms there are ^ N such interactions. The energy of the entire system may be 
written as ^ N times the ensemble average of the exchange energy 

E =1n(-2J(S 0 - Sj}) (28) 


where the angular brackets indicate ensemble averages and Hq = 0. Now, since 
2 Sq* Sj = S(S + 1) - Sj(Sj + 1) - Sq(Sq + 1) with Sq, Sp and S defined as before and 
since Sq = 1/2, then 

2( Sq • Sj) = ^S(S + 1) - + 1) - ^ 

Thus, taking the ensemble average of the exchange term 


2 (Si 


• S,> = Z -1 ^ 


y' r 

k m 


= Z” ^ ^ u>(2n, Sj) LU S(S + 1) - S 1 (S 1 



-W,/iT 
e K 


where 



u>(2n, S x ) 



k m 


-W k /iT 


results in the energy of the system having the form 



w(2n, SA 



A k e 


“W t / k T 


m 


(29) 


with all symbols as defined previously. 

As with the discussion following equation (19), the summand in the region near the 
Curie temperature may be expanded and the energy per atom £ may be obtained 
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-V* T X b k + c k |H 


A k e ^ <|1 


k m 


1 (‘■t + c k> d k 

2 (&T) 2 ^ T . 


H 


or 


iz- i y^, s 1 , yj y 


-A k /^T 


k m 



-Afc/^T 


1 ^ b k + c k^ _ d k 

2 (kr ) 2 


Hi 


For convenience in computation the energy may be rewritten in terms of four new 
functions of the temperature to second order in as 


£ = 1 ^l + l H 2 p P 1 P 2 
j‘2P 0 2 1 (p„' p 2 


( 30 ) 


where Z has been rewritten as 


z = p o + P 2 HJ 


and Pq, Pj, P 2 , and Pg are defined as 


p° ^ l yy . 


-a, At 


( 31 ) 


k m 
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-A k /iT 


(32) 


,-5> s i> 




m 




1 < b k + c k) d k 


-aJAt 


k m 


(^T) 


2 


(33) 


Si km 


~A k /AT 


(34) 


A further solution requires an explicit evaluation of these functions for the specific lattice 
of interest, that is, 2n = 5, 6, or 8. 

In addition to the temperature dependence, the energy equation (eq. (30)) contains the 
internal field explicitly. This may be eliminated by use of the self-consistency condition 
(eq. (19)) giving 



A 

C 


From these results the energy per atom has been calculated numerically as a function of 
the reduced temperature as displayed in figure 3 for the two-layer lattice, the simple 
cubic lattice, and the body-centered cubic lattice, very near the Curie temperature. It 
is evident, both from the expression for the energy (eq. (30)) and from the curves, that 
the energy is a continuous function of the reduced temperature but with a discontinuity in 


Tfl- C 


95 . 96 . 97 . 98 . 99 1.00 1.01 1.02 1.03 1.04 1.05 
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slope at the transition temperature. 

Since the specific heat of the system is simply the temperature derivative of the en- 
ergy, there is a finite discontinuity in the specific heat at the Curie temperature, as 
displayed in figure 4. The specific heat discontinuity implies that this system undergoes 
a second-order phase transition, according to the Ehrenfest definition, at the Curie tem- 
perature. The order parameter in the system is the internal field, just as in the original 
Weiss molecular field formulation. The values of the discontinuities are 0. 912, 1.211, 
and 1. 488 for the two-layered lattice, the simple cubic lattice, and the body-centered 
cubic lattice, respectively. 


Paramagnetic Susceptibility 

As is readily seen from the graph of the spontaneous magnetization for temperatures 
above the Curie temperature, the internal field vanishes in the absence of an applied 
field. When there is an applied field present, the situation is not the same. The applied 
field interacting with the spins produces an internal field and a resulting magnetization. 
This is experimentally verified by the observation that above the Curie temperature there 
is a paramagnetic susceptibility which obeys a modified Curie-Weiss law of the form 
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( 35 ) 


t - e 

where C is a constant, the Curie constant, and 9 is called the paramagnetic Curie tem- 
perature. Typically, 9 is found to be a few percent larger than T , the ferromagnetic 
Curie temperature. 

The paramagnetic susceptibility may be calculated from the BPW formalism by start- 
ing with the partition function given by equation (11), but keeping a nonzero applied field. 
For an expansion very near the Curie temperature only terms that are second order in 
the field strengths need be kept. Just as before, the average magnetic moment of the 
central atom in the cluster and the average magnetic moment of an atom in the first shell 
of nearest neighbors may be calculated from equations (15) and (16). The self-consistency 
condition requires that these be equal, giving now an equation between the applied field 
Hq and the internal field Hj 


AH + BH q = 0 (36) 

where H is the resultant internal field, H = - Hq, A is given by equation (20), and 

B has the form 


B 



b l + (1 + 2n)b k c k 
(1 + 2n)^T 



(37) 


In the calculation of the susceptibility, it is more convenient to work with the average 
magnetic moment of the entire cluster m 


m = nig + 2nm^ 


from which the magnetization m^ is 


m - m 

o 


1 + 2n 


after the self-consistency condition has been applied. 


(38) 
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( 39 ) 


k TZ 

g/3 


— =\ u(2n, Sl ) 

H 0 / ^ 



4 + b k< b k + c k>^ L 

H 0 


-Ajkl! 


The susceptibility then is the ratio of the magnetization uiq to the applied magnetic field 
Bq (or Hq/ g/3) . Thus from equations (39) and (38) 


X s 


g ^0 _ g 2 /3 2 


H f 


k TZ 


t^E w EE [ b ‘ + ^ 


+ Cfc) 


H 

Hn 


-AjkT 


(40) 


m. 


The applied field Hq and the internal field H enter only as a ratio, and this may be 
eliminated by use of the consistency equation (36) . Equation (40) may be written as 


k T 

„ 2 fi 2 
g P 




f£*(2n, Sl ) 

S, k 


II' 


m. 


•k^k + c k ) e 


■A k /kT 


(41) 


For a comparison with the Curie-Weiss law (eq. (35)), it is more convenient to cal- 
culate the reciprocal of the susceptibility. This curve then should be linear with an in- 
tercept on the T-axis that is equal to the paramagnetic Curie temperature 0. If equa- 
tion (41) is rewritten to correspond with this discussion, 


, 2 ft 2 

g P 

k T 



( 42 ) 


where Xq and x H are the corresponding summations in equation (41). Results have 
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m c 


(a) For 2n = 5. (b) For 2n = 6. 



(c) For 2n = 8. 


Figure 5. - Reciprocal of paramagnetic susceptibility as function of reduced temperature. 


been calculated for 2n = 5, 6, and 8, and are plotted in figures 5(a), (b), and (c), re- 
spectively. 

A numerical investigation of the various terms in equation (42) shows that the zero 
of the reciprocal susceptibility occurs when the factor A vanishes. But this point has 
already been found to define the ferromagnetic Curie temperature. Thus, what could be 
defined as the rigorous paramagnetic Curie temperature coincides with the ferromagnetic 
Curie temperature, at least in the cases considered here with the BPW approximation. 
Figure 5 shows, however, that another paramagnetic Curie temperature can be defined 
by extrapolating the linear portion of the curves down to the axis. The resulting values 
for the Curie temperature obtained by this scheme are 1. 36, 1. 26, and 1. 12 times the 
ferromagnetic Curie temperature for the two-layer lattice, the simple cubic lattice, and 
the body- centered cubic lattice, respectively. These extrapolated numbers do not have 
much significance due to the remaining nonlinearity in the curves for large T/T c (as 
seen in fig. 5). 


DISCUSSION 

To fulfill the original purpose for this study, it is appropriate at this point to make 
a comparison with the spin- wave results of Bloch (ref. 1), which predict a ferromagnetic 
Curie temperature for three-dimensional lattices, but none for two-dimensional systems. 
The results of the calculations presented here, applicable only to infinite lattices, seem 
to essentially agree with the results of Bloch. The linear chain and plane quadratic lat- 
tices show no ferromagnetic behavior, while the three-dimensional lattices with five, six, 

21 



and eight nearest neighbors have solutions for the Curie temperature. The Curie temper- 
ature for the two-layer thin- film model occurs at approximately two-thirds the value of 
the Curie temperature for a simple cubic lattice, both calculated with the BPW approxi- 
mation. 

The correspondence between the results of the two methods is not as complete as it 
first appears. The requirement for a lattice to have ferromagnetic properties as a result 
of the BPW method is of a more topological nature than the corresponding criterion for 
the Bloch theory. According to the BPW method, if there are no interactions between the 
first- shell atoms, any lattice that has a coordination number (number of nearest neigh- 
bors) of more than four has a Curie temperature. A calculation by Weiss (ref. 9), shows 
that lattices with interactions between the atoms in the first shell must have a coordina- 
tion number greater than or equal to eight (e. g. , a modified body-centered cubic) to pos- 
sess ferromagnetic behavior. 

As a point of contrast between the Bloch and BPW methods consider the two-layer 
model under consideration in this report, that is, a coordination number of five. It has 
been shown that such a lattice has a Curie temperature, as indicated by the nonzero solu- 
tion to equation (27). Although this represents a three-dimensional lattice, it is of inter- 
est to consider an imaginary two-dimensional lattice having the same coordination num- 
ber. The fact that such a lattice does not satisfy the symmetry requirements for a two- 
dimensional Bravais lattice is not important in this discussion. If such a plane lattice 
with coordination number five could be constructed, it would have a Curie temperature 
identical to that calculated for the two-layer lattice, that is, it would also satisfy equa- 
tion (27). 

From this example and the discussion earlier, it appears that the requirement for a 
lattice to be ferromagnetic on the BPW approximation may be stated in a concise form: 

For a lattice of a particular class (i. e. , presence or absence of interactions among the 
first- shell atoms) to possess ferromagnetic properties it is necessary and sufficient that 
the lattice coordination number be equal to or greater than some minimum value (e. g. , 

2n > 5 for noninteracting lattices). 

One point for comparison of the BPW approximation with a real ferromagnetic sys- 
tem arises from equation (30) . From this equation it is seen that the energy per atom of 
such a system is a continuous function of the temperature with a finite jump discontinuity 
in the slope at the Curie temperature. There is some evidence, on the other hand, to in- 
dicate that there should actually be a logarithmic singularity in the slope (for a discussion 
see Hill (ref. 11)). As yet, however, the experiments have not been precise enough to 
settle definitely the question of the nature of the singularity at the ferromagnetic Curie 
point. 

It is of interest to compare the results of this calculation using the BPW approxima- 
tion with the "constant- coupling" approximation of P. W. Kastelijn (ref. 12). This 
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method assumes that the Heisenberg Hamiltonian for the entire system may be written in 
terms of effective scalar two-body interaction Hamiltonians. From these two-body terms 
the energy of the system may be calculated by means of standard thermodynamic tech- 
niques. 

The interest in this method arises from the fact that for a cluster that contains no 
interactions between first- shell atoms, the BPW method represents, in a sense, a sys- 
tem of two-body interactions, just as in the Kastelijn approximation. This correspon- 
dence is not complete, since the BPW method includes long-range correlations inherent 
in the assumption about the internal field, not included in the constant- coupling method. 
For the two-layer lattice, however, first- shell interactions are absent, and the constant- 
coupling method gives a useful cross-check for the calculations done previously. 

The expressions derived by Kastelijn for the Curie temperature T and the discon- 

v 

tinuity in the specific heat AC are 

At = 2jYln 2n -j 1 
c V 2n - 4/ 

AC = — 4n2 (^- l ) 2 .. (i n _2n_fj 

32 (2n - l)(2n - 2) \ 2n - 4/ 

where 2n represents the number of nearest neighbors for the lattice, as used before. 
Evaluating these quantities for 2n = 5 results in 

T = 1. 243 — 

C k 

AC = 0. 506 k 

which compare fairly well with our value of 1. 278 ^ for the Curie temperature, but not 
with the specific heat discontinuity (0. 912 A) obtained from the BPW approximation. 

There are two points to be noted in making a comparison between the BPW and the 
constant- coupling approximations. First, the condition that the lattice must have 2n > 5 
for ferromagnetic behavior to exist is the same in each method. Secondly, in the 
constant- coupling approximation there is no second solution corresponding to the anti- 
Curie point as appears in the BPW calculations. 

CONCLUDING REMARKS 

Within the limitations of the Bethe-Peierls-Weiss approximation the two-layer lattice 
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model for a thin film shows behavior that is qualitatively consistent with the results of 
spin-wave theory and also with the properties of noninteracting bulk lattices as calculated 
by the BPW method. The Curie temperature is consistent with that calculated for the 
simple cubic and body- centered cubic, and the second-order nature of the phase transi- 
tion appears as a finite jump discontinuity in the specific heat. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, June 1, 1965. 
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APPENDIX - CLUSTER HAMILTONIAN 


This appendix deals with the problem of calculating the matrix elements of the 
Hamiltonian (eq. (3)) in a representation in which the exchange energy is diagonal. The 
discussion essentially follows that of Condon and Shortley (ref. 13). The Hamiltonian 
may be written in terms of Hq and Hj as discussed in the text 

^ = -2JV S, - H 0 S 0z - Hl S lz (Al) 

where the vector sum of S Q and Sj is defined by 

S = S Q + S x (A2a) 


and 


S z = S 0z 


+ S 


lz 


(A2b) 


Then, by means of equations (4) and (A2b) dt may be rewritten in a more useful form 


<# = -J 




S Z H 1 + S 0z< H l - H 0> 


(A3) 


The representation that diagonalizes the exchange interaction is the representation 
in which S , S z , Sq, and Sj are diagonal. The eigenvalues of the exchange energy are 
given by equation (6) and the eigenvalues of S„ are just the values of m. All that re- 
mains is to determine the matrix elements of Sq z in this representation. 

To accomplish this, consider the properties of a class of operators known as 
"class T M operators. These operators are defined by their commutation relations with 
any angular momentum operator J. These relations are 


o 

il 

J 

[j y> T x ]=-iT z 

[J I( TJ- 1T y ' 


[J X , = iT z 

f J y’ T y^ = 0 

[ J z’ T yJ = - 1T x 

► (A4) 

r— i 

X* 

£ 

It 

[^y> T z ] = iT x 

O 

II 

N 


and 





[j,T 1 -T 2 ] = o 


(A5) 
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where square brackets symbolize commutator brackets and K = 1. Any vector operator 
that transforms in the same manner as the position vector under a proper coordinate ro- 
tation is a class T operator. In particular, an angular momentum operator, such as 
Sq, belongs to this class of operators. 

To calculate the matrix elements of Sq z , plan on taking the matrix elements of one 
component of a vector operator T. For this purpose, consider 

[j 2 ,T] = J- [J,T] - [T, J] • J 


= -i(J X T - T X J) 


= -2i(J X T - iT) 


which follow from the commutation relations. Now consider the following commutator: 

[J 2 , [J 2 , T]] = -2i[J 2 , J x T - if] 

= -2i[J 2 , J x T] - 2[J 2 , T] 

= 2 ( J 2 T 2 + TJ 2 ) - 4J(J • T) 

But there is also another expression for [J , [J , T]] 

[j 2 , [J 2 ,T]] = J 4 T - 2J 2 TJ 2 + TJ 4 


Therefore, 


J 4 T - 2J 2 TJ 2 + TJ 4 = 2(J 2 T + TJ 2 ) - 4J(J- T) (A6) 

From this identity the matrix elements of T can be calculated since J is a Her- 
mitian operator. For example, from the left side of equation (A6), for the matrix ele- 
ments between the states (jm| and Ij’m'), as 

(jm|j 4 T - 2J 2 TJ 2 + TJ 4 1 j'm’) = (jm | J 4 T [ j'm*) - 2(jm | J 2 TJ 2 Ij'm') 

+ ( jm |TJ 4 | j'm') = [j(j + 1) - j'(j* + l)] 2 (jm|f |j’m’) (A7) 

since a system of units is used in which K = 1. Similarly for the ’right side of equation 
(A6), 
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(jm|2(J 2 T + TJ 2 ) - 4J(J • T) | j’m’) = 2[j(j + 1) + j'(j* + 1)] (jm|Tj j'm*) 

-4{ jmjj(J- T) Jj'm') (A8) 


But the right side of equation (A8) may be rewritten as 


T) | j’m') 


<jm|j|j>*m’*)<j"m"|j- Tjj'm') 

j'*m" 


But J • T is diagonal with respect to j and m, which follows from equation (A5); thus 

(jm)j(J- T) j j'm*) = (jm|j)jm’><jm*|j- T|jm*> 

Thus, 

|[j(j + 1) ~ j'(j* + l)] 2 - 2[j(j + 1) + j'(j' + l)]|{jm[T|j , m*> 

= -4(jm|j|jm')(jm'|j- T|j'm*) (A9) 

This is the expression determining the matrix elements of T, and in particular, those 
of T z . 

Note in equation (A9) that, if j’ = j ± 1, the square bracket on the left side vanishes. 
Looking only for the expression when j = j’ gives 


j(j + 1){ jm | T | jm’) = (jm | j| jm') (jm'| J • T|jm’} 


(A10) 


Since J z is diagonal in m, it follows for the z- component that 

j(j + l)(jm|T z |jm) = (jm | J z ] jm) ( jm | J • T|jm) (All) 


This expression determines the diagonal elements of S 0z in this representation. 

Since equation (A9) is indeterminate for the off-diagonal elements, another expres- 
sion must be derived that is valid when j* = j ± 1. To do this, consider the ladder oper- 
ators T + and T_ defined on the set of class T operators 


T = T + iT 
+ x y 


T = T - iT 
x y 


(A12) 
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Similar ladder operators may be defined for the J operators, as 


J + " J x + iJ y 


J -= J x- iJ y- 


(A13) 


To express T in terms of these operators, consider the following commutator: 
z 


[J + ,TJ=[J x + U y ,T x -lT] 


= [J., T_] + 2i[Jy, TJ 


or, from equation (A4), 


T z =i[J + ,T_] (A 14) 

2 

Thus, from this expression and a knowledge of the matrix elements of the ladder oper- 
ators the matrix elements of T z may be calculated. For example, 


<jm|T z |j'm)=l ({jm|j + |j”m»«)<j'»m»*|T_|j*m>-<jm|T_|j**m»*>(j»m*»|j + |j»m)) 


j"m" 


(A15) 


The ladder operators act only upon the m values and have nonzero matrix elements only 
between two states whose m-values differ by unity. 

Consider each case in equation (A15) separately (see refs. 13 to 15) 

(1) For j* = j - 1 

2<jm|T z [j - 1, m) = <jm|j + |j,m - l)(j,m - l[T_|j - 1, m) 

- <jm|T_ |j - l,m + l)(j - l,m + l|j + |j - 1, m) 

= - 1) V(j - m + l)(j - m) V(j + m)(j - m + 1) 


= 2 TjZ 


m 2 (j:T:j - 1) 


- (j:T:j - 1) V(j - m)(j - m - 1) V (j + m)(j - m + 1) 

(A16) 
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(2) For j* = j + 1 


2<jm|T z |j + l,m> = <jm|j + |j,m - l><j,m - 1 1 T_ [ j + l,m) 

- (jm |T_ | j + l,m + l)(j + 1, m + l|j + [j + 1, m) 

= 2f(j + l) 2 - m 2 (j:T:j + 1) (A17) 

where the quantities (j:T: j - 1) and (j:Ti j + 1) are in the notation of Condon and 
Shortley (ref. 13) and need to be defined. Starting with the commutation relation for J_ 
and T_, that is, [J_, T_] = 0, and calculating the matrix elements that satisfy the selec- 
tion rules on m result in 

<j,m - l|j_|jm)(jm|T_|j , ,m + 1) = (j,m - 1 |T_ |j'm)(j , m| J_ | j», m + 1) (A18) 

for j* - j = 0, ±1. 

But the matrix elements of the J_ operator are well-known and yield 

V(j +m)(j - m + 1) (jm |t_ | j’, m + 1) = V(j* + m + l)(j* - m) (j,m - l|T_|j*m) 


Thus 

<Jm|T_|j',m + 1) ( j> m - 1 [t_ | j', m) 

f j’ - j = 0, ±1 (A19) 

Vo* + m + 1) (j * - m) Vo" + m) (j - m + 1) 


This relation between these ratios is independent of m. 

This equality evaluated for j* = j ± 1 is used to define the terms (j:T: j - 1) 
<j:T:j + 1) 

( jm |T_ | j - 1, m + 1) <j,m - l|T_|j - l,m> 

= — ■■■ - ■ = <j:T:j - 1) 

V(j - m - 1) (j - m) V(j - m) (j - m + 1) 

► 

(jm |t_ | j + 1, m + 1) (j,m - 1 1 T_ | j + 1, m) 
-—^^=^=2 = — ' = "(i:T:j + 1) 

V(j + m + 2)(j + m + 1) \(j + m + 1)0 + m) -J 


and 


(A20) 


These elements are related to the set of irreducible elements of a class T operator. 
The evaluation of these will be reserved until later. 
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Equations (All), (A16), and (A17) determine the matrix elements of Sq z explicitly, 

once a change of notation is made. Before this is done, it would be useful to consider the 

general form of the operator Sq z . This matrix will have 2 + 1 = 2(S^ + 1) rows 

and columns. The only nonzero elements will lie upon the principal diagonal and along 

the diagonals one step removed from the principal diagonal. This follows from the fact 

that the only nonzero matrix elements for T are for j' - j = ±1. Each row or column is 

indexed by a value of Sj, or, equivalently by m. Calculating the matrix elements of 

S~ as given in the text results in effectively picking out a square 2 by 2 matrix with 
uz i i 

j = ®1 + 2’ ®i “ "5 for a given value of m. Thus to calculate the matrix elements 
from equation (All), (A16), and (A17), the association is made that the j-values for the 
two states are Sj + and Sj - associating with T z the operator S 0z> which is a 
class T operator. 

First, for the diagonal elements of Sq z , (jm|j - T|jm) is needed for T — Sq and 
J — S. Consider Sj = S - S^; then 

S 2 = S 2 + - 2S q - S 

or, 

S 0 S = i(s 2 + s2-s2) 

and then, 

4> m iv i ii s i 4' m ) 4 [( s i 4)( Si + !H ■ Si(Si + J) ] 4 ( Si + 2 



and 


Thus, 



(A21) 
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m 

2Sj + 1 


(A22) 


For the off-diagonal matrix elements the irreducible terms appearing in equations 
(A16) and (A17) must be evaluated. When Sq is the operator T involved, the evaluation 
of these terms may be derived, but the calculation is prohibitively long to be shown here. 
The details are given in Condon and Shortley (ref. 13, chapter 3); the results are merely 
quoted here 




(S - S 0 + SjMS + S 0 - S 1 )(S + S Q + S 1+ 1)(S 0 + S x + 1 - S) 


4S^(2S - 1)(2S + 1) 


1/2 


for 


8=8, +- 

1 2 


and 




(S + 1 - S 0 + S^S + 1 + S 0 - S 1 )(S 0 + S x - S)(S 0 + S x + S + 2) 
{ 4 (S + 1) 2 [4(S + l) 2 - l]} 


1 1/2 


for S = S, - — . 

1 2 

Evaluating these gives 




•S 0 :Si 


i)— i— 

2/ 2Sj + 1 


(A23) 


and 




1 

2S 1 + 1 


(A24) 


Substituting these results into equations (A16) and (A17), the matrix elements of Sq z 

are 
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^ S 1 +^ m |S 0z |Si ^ - 2 ’ m l S 0zl S l + 2 ’ m ) 2 


1 - 


4m" 


(2S 1 + 1)' 


1/2 


Thus, from equations (A21), (A22), and (A25), the matrix form of one of the 2 by 2 

1 / 2 ' 


block matrices of Sq z is 


f - s - 

/ 2Sj + 1 


®0z 


1 - 


4m 


L (2S 1 + 1) J 


(A26) 


V 


4m 


(2S 1 + if] 


1/2 


m 


2S X + 1 


Accordingly, the matrix of the Hamiltonian is 


JSf = 


f J<Sj + 1 ) - -H - (Hj - Hq) 


4m" 


(2Sj + 1)" 


1/2 


(H 


i-% \ 


V 


4m" 


( 2 Sj + if] 


1/2 


(Hi - H q ) 


-JS^ - mHj + 


m 


2Sj + 1 


(Hi - H q ) 


(A27) 


as given in the text. These are the results given by Weiss (ref. 9). 
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